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Abstract. Electrons released from clusters through strong Xray pulses show broad kinetic- 
energy spectra, extending from the atomic excess energy down to the threshold, where usually a 
strong peak appears. These low-energy electrons are normally attributed to evaporation from the 
nano-plasma formed in the highly-charged clusters. Here, it is shown that also directly emitted 
photo electrons generate a pronounced spectral feature close to threshold. Furthermore, we give 
an analytical approximation for the direct photo-electron spectrum. 
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1. Introduction 

Recently, there has been increasing interest in slow electrons from photo-driven 
processes. While identified in strong-held ionization of atoms [1,2] as well as molecules 
[3], the mechanisms behind the production of slow electrons are very different for 
atoms in linearly polarized pulses [4] and molecules in elliptically polarized pulses [5], 
respectively. Crucial in both cases is the (single) electron dynamics in the combined 
potential of the ion left behind and the driving laser held. 

Slow electrons can also emerge from soft and even hard Xray pulses. At a hrst 
glance this is surprising, since the electronic excess energy E* (which is the photon 
energy reduced by the binding energy) is typically large, say a few hundred eV up to 
few keV, depending on the photon energy. Under such circumstances, the low-energy 
electrons can occur through non-adiabatic effects in very short pulses, when the pulse 
length becomes comparable with the orbital period of the bound orbital which is photo- 
ionized [6,7]. 

While this effect is again essentially a single-electron phenomenon, another very 
common mechanism to produce slow electrons in intense Xray pulses requires although 
not collective, yet multiple ionization: Thereby, a complex of ions (either clusters or 
big molecules) staying behind forms a large background charge [8-15], which reduces 
the excess energy E*. Hence, these kind of slow electrons can only emerge from large 
systems, which allow for high charging. In fact, the background charge may be so 
large that electrons are being trapped even for photons in the keV-range [10]. The 
trapping leads to the formation of a so-called nano-plasma, which thermalizes quickly 
and consequently evaporates (slow) electrons. Typically, the yield of the slow electrons 
shows an exponential decrease with an energy-scaling constant related to the plasma 
temperature according to common sense. However, this relation is tricky for two reasons. 
Firstly, due to the continuous excitation of electrons into the plasma its state may change 
considerably during the Xray pulse violating the quasi-stationarity which is necessary to 
assign a temperature to the electron spectrum. Secondly, the photo-ionization process 
itself gives rise to directly ejected slow electrons. 

This is illustrated in Fig. 1 with the electron spectrum resulting from the 
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Figure 1. Electron spectra calculated for a Coulomb 
complex [13] of radius R = 10 with 10 3 electrons for an 
excess energy of E* = 50 and a pulse duration of T = 32. 
Details of the numerical approach are given in section 
4. The full spectrum is shown as darkgray-shaded area, 
the one for direct electrons only, i. e. excluding plasma- 
electrons, as red solid line. 
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illumination of a generic (spherical) cluster with a short pulse with T = 32 duration 
and excess energy of E* = 50. The direct electrons (red curve) show a clear peak at 
low energies revealing that the slow electrons do not only result from the evaporation 
of the nano-plasma. We define the direct electrons as those electrons which have a 
positive energy p* 2 /2 + W( r*) > 0 just after the absorption of a photon at time t* 
at position r* = r (t*) with momentum p* = p(f*). The potential energy W involves 
both the attraction from the ionic background, defined below in section 4, and the 
electron-electron repulsion. The initial momentum p* is fixed by the excess energy E*. 

In order to understand the peak in the (numerically obtained) direct-electron yield 
in Fig. 1, we will provide in section 2 an analytical derivation of the direct-electron 
spectrum under the premise that the direct electrons leave the cluster sequentially and 
(indirect) plasma electrons remain in the cluster. Thereby, the origin of the slow direct 
electrons will become clear. With a surprisingly simple approximation, suggested by the 
form of the direct-electron spectrum, we can give a fully analytical formula (section 3). 
It is compared in section 4 to the numerical direct spectrum, revealing how the indirect 
plasma electrons influence the direct electrons. 

2. The direct photo-electron spectrum and the origin of slow direct 
electrons 

We assume here for simplicity that the system is spherical with a radius R throughout 
the ionization process. The light pulse leads to random single-ionization events of atoms 
within the cluster, where we choose the intensity such that the system is far from 
saturation of complete single ionization and the occurrence of any multiple ionization of 
cluster atoms. More explicitly, if the cluster contains N atoms and the pulse leads to Q 
ionization events, then in the end N—Q « N atoms of the cluster remain neutral. For 
the case of sequential ionization the photo-electron spectrum follows from integrating 
the spectra P q for an instantaneous charge q ranging from q = 0 (for the initially neutral 
cluster) to q = Q (the highest possible charge state) 

P(E)= f AqP,(E). (1) 

Jo 

The highest charge Q is reached when the cluster potential V q (r) = [q/2R][3 — r 2 /R 2 ] is 
so deep that absorption of a single photon (with excess energy E*) is not sufficient to 
overcome the threshold*. This occurs if Vq(R) = E* which implies Q = E*R. 

If the cluster potential V q (r) is still shallow enough for all electrons in the cluster to 
escape by absorbing just one photon, the electron spectrum of a (/-fold charged spherical 
cluster is given by [17] 

p ,(E) = -^ J drr 2 S(E-E,(r)) (2a) 

f Note that higher charge states can be reached when electrons are excited below threshold and the 
nano-plasma, formed in the process, evaporates [8,11,16]. 
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Figure 2. Sketch of formation of the photo-electron spectrum (a) Final spectrum as obtained by 
numerical integration (blue dashed line) of Eq. (1) with Eqs. (3) and (4) and from the analytical 
approximation (red solid line) according to Eq. (9); (b) The contribution from a particular charge: 
The shaded area shows which charges q contribute to which energy E either according to P q 
(light-gray) or P* rap (dark-gray), respectively. Additionally, there are three explicit examples 
with q = 2Q/5, 3Q/5, 4Q/5 for these distributions according to Eqs. (3,4) and (8), respectively. 
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the final energy of an electron released at a distance r from the centre through the 
absorption of a photon. We obtain from Eq. (2) 


P q (E) 


3 I E*-E 

q/R\ 2 q/R 

with E min (q) = E* 


for E min (q) < E < E max (q) 
3q/2R and E max (q ) = E* - q/R 


(3a) 

(3b) 


and P q (E ) = 0 elsewhere. Here, E min (q ) = E q ( 0) is the energy from an electron released 
at the centre (r = 0), while an electron from the surface will appear at E max (q) = E q (R ), 
see Eq. (2b). The two lower blue dashed lines in Fig. 2b show as examples P q {E) for 
q=2Q/5 and q=3Q/5, respectively. 

Expression (3) has to be modified when the cluster potential becomes so deep 
that electrons — firstly those released at the centre — are trapped after single-photon 
absorption. This occurs at g trap = E*2R/3 = 2Q/3. Since these electrons do not 
increase the cluster charge, further charging is due to electrons located closer and closer 
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to the surface. Hence, for charges q tIap < q < Q the instantaneous spectra become 


P,(E) -> P^(E) =- 1 „ -P,(g) for 0 <E<E m „(g) (4) 

^{^m\ 


and Pq rap (E) = 0 elsewhere. The prefactor in (4) normalizes the distribution P* rap for 
any q just as P q above is normalised, i.e., 

pE* pE* 

/ d EP q {E)= / dEP^ p (E) = 1. (5) 

Jo Jo 

One example for P* rap (P) is shown in Fig. 2b with the upper blue dashed line 
corresponding to q—4Q/5. One can also see from the shaded area in Fig. 2b that the 
restriction of electron energies to the interval E mm (q) < E(q) < E max (q) implies for the 
integral (1) a restriction to charges in the interval q m i n (E ) < q(E) < q ma x(E ) with 




-[E*-E]R and 
o 


Qm&xi.E') 


[E*-E\R, 


( 6 ) 


which follows directly from Eq. (3b). The abundance for a particular energy E finally 
reads 


'Qr nax 


P(E )=/ d qP q (E) 

for E */3 < E < E* 

(7a) 

** Qmin 

pq trap Nmax 



P(E)= d qP,(E)+ d qP'^(E) 

for 0 < E < E*/ 3. 

(7b) 


Q'min 'J (/trap 


Equation (7a) can be solved analytically and gives the energy-independent value 
P[E ) = 3R[VS ln(2+V3)-2] corresponding to a plateau [13]. Equation (7b), on 
the other hand, does not allow for a compact analytically solution. Therefore, we 
provide with the blue dashed line in Fig. 2a the numerically integrated spectrum. One 
clearly sees an accumulation towards lower energies with a divergence at E = 0. With 
(/ min and (/trap finite, this is due to the second term in (7b) and may be interpreted as 
follows: Electrons with energies E > E* /3 can escape from anywhere in the cluster 
for any (accessible) charge state, as described by Eq. (7a). For electrons with energies 
E < E* /3 this is limited to clusters charged less than g trap = 2Q/3, cf. 1st integral in 
(7b). For clusters charged higher than (/trap direct electrons come from the outer regions 
of the cluster with ever decreasing energy as the cluster charge grows beyond (/trap- This 
part of the spectrum is described by the 2nd integral in (7b). 


3. Analytical approximation for the direct photo-electron spectrum 

Interestingly, the exact shape of P q (E) is not important for the final spectrum. One 
may choose any form for P q . As long as the shape for various values q can be obtained 
by a simple scaling the final spectrum is a plateau [17]. In order to obtain an analytical 
expression for all energies we approximate P q (E) with the simplest form possible, namely 
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a constant spectrum between E min and E max . The (/-dependence is then introduced via 
the normalization (5). With E min (q) and E max (q ) given in (3b) this leads to 


P,(E) 

P,(E) 


2 R 


^max -^min Q 


E n 


E* — q/R 


for 

0 <q< i/trap 

(8a.) 

for 

(/trap Q — Q ■ 

(8b) 


This distributions are shown in Fig. 2b as red solid lines. They allow for an integration 
of Eq. (7) 


P(E) = 2 R ln(3/2) for E */3 < E < E* (9a) 

7 7 

P(E) = 2R In ( WZE ) + R In (^) for 0 < E < E'/3 (9b) 

which is shown as red solid line in Fig. 2a. ft reproduces the spectrum obtained with 
the blade-shaped instantaneous spectra (blue dashed line in Fig. 2a) extremely well. 
This applies to both, the absolute values of the plateau at large E, and the divergent 
behaviour around E— 0. Note, that the latter is indeed due to the second term in (9b) 
which represents the analytical approximation of the integral over P 4rap in (7b). 


4. Comparison to numerical results from Coulomb complexes 

The analytical expressions have been derived under the assumptions that photo¬ 
ionization occurs sequentially and that electrons excited to states below threshold remain 
trapped. In the following we will assess if and when these assumptions are justified by 
comparing the results from (9) to those of molecular dynamics calculations without 
those assumptions in the framework of so-called photo-activated Coulomb complexes 
[13]. This is a simple model, where electrons are treated as classical particles and ions 
form a spherical jellium, describing the attractive potential of the charged cluster. This 
potential as well as the electron-electron interaction is essential for understanding the 
formation of the broad electron spectra. Once activated (i. e. released with a given 
energy determined by the excess energy E*) the electrons are propagated according to 
Newton’s equations with forces resulting from the jellium potential and the electron- 
electron interaction. Thus, in contrast to the description of the previous section, here 
correlations (collisions) of the electrons are fully taken into account. The system is 
propagated sufficiently long (up to times t = 10 4 for the results presented) before spectra 
are calculated. These are obtained by folding the final (kinetic) energies Ej of the 
electrons with a Gaussian 

P(E) = ^ exp (-[Ei-EY/SE*) (10) 

3 

of width 5E= 1. 

Figures 1 and 3 show such spectra (obtained by averaging over 100 realizations) 
for a Coulomb complex of radius R = 10 with 10 3 electrons and E* = 50. The photo 
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Figure 3. Electron spectra (red solid line) for the direct photo-electrons only, i. e. excluding 
plasma-electrons, as obtained from Coulomb complexes with 10 3 electrons for an excess energy 
E* = 50 and various pulse durations T. They should be compared to the analytical expressions 
(9), which is shown by the gray-shaded areas. 


activation rate is proportional to exp (—t 2 /T 2 ). One clearly sees a broad spectrum with 
a large peak at E ps 0, a plateau at E < E* and a cutoff at E = E*. These features 
have been observed [18] and discussed [11-14] before, interpreting the high-energy part 
(plateau) as a consequence of the direct photo electrons and the low-energy part with 
its peak towards threshold as a consequence of the evaporation from the transient nano- 
plasma. 

However, as already mentioned, Fig. 1 reveals that also direct electrons, defined as 
those electrons which have initially enough energy to escape from the cluster potential, 
contribute to the slow-electron peak. Their contribution to the low energy spectrum 
is even larger in the analytical estimate considering only sequentially emitted electrons 
(see Fig. 3, gray-filled area) than from the numerically obtained direct electrons (red 
curve). The reason is that we do not take into account that initially trapped plasma 
electrons do eventually leave at a certain rate, dictated by the plasma temperature. If 
this rate is faster than the photo ionisation rate, direct photo electrons see an increased 
background charge reducing their yield at low energies since they get trapped. This 
effect should be least important for very short pulses when the direct electrons leave 
before plasma evaporation becomes important. However, for very short pulses, the 
second assumption made for the analytical direct electron spectrum is violated, namely, 
the sequential ionization: The photo-ionization rate is so large that the direct electrons 
interact and exchange energy before leaving the cluster. This indicates the onset of 
massively parallel ionization [19], which is accompanied by high-energy tails at E>E* 
in the spectrum. Indeed, the red curves from the numerical calculation in Fig. 3 show 
these tails in contrast to the sharp cutoff of the analytical spectrum at E = E*. 












5. Summary 


Comparing fully numerical spectra to those from photo electrons only, we have shown 
that the low-energy peak observed in the photo-electron spectrum of multiple ionization 
of clusters in strong Xray pulses is not only generated by initially trapped plasma 
electrons but also by photo electrons directly escaping. An understanding of the origin 
of slow direct electrons has been made possible by the formulation of the spectrum for the 
direct electrons alone down to threshold, including a fully analytical approximation — 
always under the assumption that the electrons leave the cluster sequentially. In 
the future it would be interesting to disentangle direct photo-electron dynamics from 
plasma-electron dynamics experimentally. This could be done by by exploiting the fact 
that angular distributions are different for photo-electrons (depending on the shape of 
the orbitals being ionized) and plasma-electrons (expected to be isotropic) or by using 
streaking techniques [20]. 
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